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The inner 0.1 AU around accreting T Tauri stars hold clues to many physical processes 
that characterize the early evolution of solar-type stars. The accretion-ejection connection 
takes place at least in part in this compact magnetized region around the central star, with 
the inner disk edge interacting with the star's magnetosphere thus leading simultaneously to 
magnetically channeled accretion flows and to high velocity winds and outflows. The magnetic 
star-disk interaction is thought to have strong implications for the angular momentum evolution 
of the central system, the inner structure of the disk, and possibly for halting the migration of 
young planets close to the stellar surface. We review here the current status of magnetic field 
measurements in T Tauri stars, the recent modeling efforts of the magnetospheric accretion 
process, including both radiative transfer and multi-D numerical simulations, and summarize 
current evidence supporting the concept of magnetically-channeled accretion in young stars. 
We also discuss the limits of the models and highlight observational results which suggest that 
the star-disk interaction is a highly dynamical and time variable process in young stars. 



1. THE MAGNETIC ACCRETION PARADIGM 

T Tauri stars are low-mass stars with an age of a few 
million years, still contracting down their Hayashi tracks to- 
wards the main sequence. Many of them, the so-called clas- 
sical T Tauri stars (CTTSs), show signs of accretion from a 
circumstellar disk (see, e.g., Menard and Bertout, 1999 for 
a review). Understanding the accretion process in T Tauri 
stars is one of the major challenges in the study of pre-main 
sequence evolution. Indeed, accretion has a significant and 
long lasting impact on the evolution of low mass stars by 
providing both mass and angular momentum. The evolution 
and ultimate fate of circumstellar accretion disks have also 
become increasingly important issues since the discovery 
of extrasolar planets and planetary systems with unexpected 
properties. Deriving the properties of young stellar systems, 
of their associated disks and outflows is therefore an impor- 
tant step towards the establishment of plausible scenarios 
for star and planet formation. 

The general paradigm of magnetically controlled accre- 
tion onto a compact object is used to explain many of the 
most fascinating objects in the Universe. This model is a 
seminal feature of low mass star formation, but it is also en- 
countered in theories explaining accretion onto white dwarf 
stars (the AM Her stars, e.g., Warner, 2004), accretion 



onto pulsars (the pulsating X-ray sources , e.g., Ghosh and 
Lamb, 1979a), and accretion onto black holes at the cen- 
ter of AGNs and microquasars (Koide et ah, 1999). Strong 
surface magnetic fields have long been suspected to exist 
in TTSs based on their powerful X-ray and centrimetric ra- 
dio emissions (Montmerle et al, 1983; Andre, 1987). Sur- 
face fields of order of 1-3 kG have recently been derived 
from Zeeman broadening measurements of CTTS photo- 
spheric lines (Johns Krull et al, 1999a, 2001; Guenther 
et al, 1999) and from the detection of electron cyclotron 
maser emission {Smith et al, 2003). These strong stellar 
magnetic fields are believed to significantly alter the accre- 
tion flow in the circumstellar disk close to the central star. 

Based on models originally developed for magnetized 
compact objects in X-ray pulsars (Ghosh and Lamb, 1979a) 
and assuming that T Tauri magnetospheres are predomi- 
nantly dipolar on the large scale, Camenzind (1990) and 
Konigl (1991) showed that the inner accretion disk is ex- 
pected to be truncated by the magnetosphere at a distance 
of a few stellar radii above the stellar surface for typical 
mass accretion rates of 10~ 9 to 10~ 7 M Q yr _1 in the disk 
(Basri and Bertout, 1989; Hartigan et al, 1995; Gullbring 
et al, 1998). Disk material is then channeled from the disk 
inner edge onto the star along the magnetic field lines, thus 
giving rise to magnetospheric accretion columns. As the 
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free falling material in the funnel flow eventually hits the 
stellar surface, accretion shocks develop near the magnetic 
poles. The basic concept of magnetospheric accretion in T 
Tauri stars is illustrated in Figure 1 . 

The successes and limits of current magnetospheric ac- 
cretion models in accounting for the observed properties of 
classical T Tauri systems are reviewed in the next sections. 
Sect. 2 summarizes the current status of magnetic field mea- 
surements in young stars, Sect. 3 provides an account of 
current radiative transfer models developed to reproduce the 
observed line profiles thought to form at least in part in ac- 
cretion funnel flows, Sect. 4 reviews current observational 
evidence for a highly dynamical magnetospheric accretion 
process in CTTSs, and Sect. 5 describes the most recent 2D 
and 3D numerical simulations of time dependent star-disk 
magnetic interaction. 




Fig. 1. — A sketch of the basic concept of magnetospheric 
accretion in T Tauri stars (from Camenzind, 1990). 



2. MAGNETIC FIELD MEASUREMENTS 

2.1. Theoretical Expectations for T Tauri Magnetic 
Fields 

While the interaction of a stellar magnetic field with an 
accretion disk is potentially very complicated (e.g., Ghosh 
and Lamb, 1979a,b), we present here some results from the 
leading treatments applied to young stars. 

The theoretical idea behind magnetospheric accretion is 
that the ram pressure of the accreting material (P ram = 
0.5/w 2 ) will at some point be offset by the magnetic pres- 
sure (Pb = B 2 /8ir) for a sufficiently strong stellar field. 
Where these two pressures are equal, if the accreting ma- 
terial is sufficiently ionized, its motion will start to be con- 
trolled by the stellar field. This point is usually referred to 
as the truncation radius (Rt)- If we consider the case of 
spherical accretion, the magnetic field becomes 



If we then assume a dipolar stellar magnetic field where 
B = £?*(i?*/r) 3 and set the velocity of the accreting ma- 
terial equal to the free-fall speed, the radius at which the 
magnetic field pressure balances the ram pressure of the ac- 
creting material is 
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where B 3 is the stellar field strength in kG, M_ 8 is the mass 
accretion rate in units of 1CP 8 M Q yr -1 , Mo. 5 is the stellar 
mass in units of 0.5 M Q , and R2 is the stellar radius in units 
of 2 R . Then, for B* = 1 kG and typical CTTS properties 
(M, = 0.5 M , R* = 2 R , and M = 10~ 8 M y r - 1 ), 
the truncation radius is about 7 stellar radii. 

In the case of disk accretion, the coefficient above is 
changed, but the scaling with the stellar and accretion pa- 
rameters remains the same. In accretion disks around young 
stars, the radial motion due to accretion is relatively low 
while the Keplerian velocity due to the orbital motion is 
only a factor of 2 1 / 2 lower than the free-fall velocity. The 
low radial velocity of the disk means that the disk densities 
are much higher than in the spherical case, so that the disk 
ram pressure is higher than the ram pressure due to spheri- 
cal free-fall accretion. As a result, the truncation radius will 
move closer to the star. In this regard, equation 2.2 gives 
an upper limit for the truncation radius. As we will dis- 
cuss below, this may be problematic when we consider the 
current observations of stellar magnetic fields. In the case 
of disk accretion, another important point in the disk is the 
corotation radius, Rco, where the Keplerian angular veloc- 
ity is equal to the stellar angular velocity. Stellar field lines 
which couple to the disk outside of Rco will act to slow 
the rotation of the star down, while field lines which couple 
to the disk inside Rco will act to spin the star up. Thus, 
the value of Rt relative to Rco is an important quantity 
in determining whether the star speeds up or slows down 
its rotation. For accretion onto the star to proceed, we have 
the relation Rt < Rco ■ This follows from the idea that at 
the truncation radius and interior to that, the disk material 
will be locked to the stellar field lines and will move at the 
same angular velocity as the star. Outside Rco the stellar 
angular velocity is greater than the Keplerian velocity, so 
that any material there which becomes locked to the stellar 
field will experience a centrifugal force that tries to fling the 
material away from the star. Only inside Rco will the net 
force allow the material to accrete onto the star. 

Traditional magnetospheric accretion theories as applied 
to stars (young stellar objects, white dwarfs, and pulsars) 
suggest that the rotation rate of the central star will be set 
by the Keplerian rotation rate in the disk near the point 
where the disk is truncated by the stellar magnetic field 
when the system is in equilibrium. Hence these theories 
are often referred to as disk locking theories. For CTTSs, 
we have a unique opportunity to test these theories since 
all the variables of the problem (stellar mass, radius, rota- 
tion rate, magnetic field, and disk accretion rate) are mea- 
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sureable in principle (see Johns-Krull and Gafford, 2002). 
Under the assumption that an equilibrium situation exists, 
Konigl (1991), Cameron and Campbell (1993), and Shu et 
al. (1994) have all analytically examined the interaction 
between a dipolar stellar magnetic field (aligned with the 
stellar rotation axis) and the surrounding accretion disk. As 
detailed in Johns-Krull et al. (1999b), one can solve for the 
surface magnetic field strength on a CTTS implied by each 
of these theories given the stellar mass, radius, rotation pe- 
riod, and accretion rate. For the work of Konigl (1991), the 
resulting equation is: 

/ e \ V6 / 3 \ -7/4 / M« \ 5/6 

In the work of Cameron and Campbell (1993) the equation 
for the stellar field is: 

/ B \ -3 / P \ 29/24 

*(fj) (rt) kG - < 2 ' 4 > 

Finally, from Shu et al. (1994), the resulting equation is: 

( a x \ -7/4 / M* \ 5/6 / M \ 1/2 

S *- 3 ' 38 1 0923 J ImqJ llO-7M yr-J X 

All these equations contain uncertain scaling parameters 
(e,3, 7, a x ) which characterize the efficiency with which 
the stellar field couples to the disk or the level of vertical 
shear in the disk. Each study presents a best estimate for 
these parameters allowing the stellar field to be estimated 
(Table 1). Observations of magnetic fields on CTTSs can 
then serve as a test of these models. 

To predict magnetic field strengths for specific CTTSs, 
we need observational estimates for certain system param- 
eters. We adopt rotation periods from Bouvier et al. (1993, 
1995) and stellar masses, radii, and mass accretion rates 
from Gullbring et al. (1998). Predictions for each analytic 
study are presented in Table 1 . Note, these field strengths 
are the equatorial values. The field at the pole will be twice 
these values and the average over the star will depend on 
the exact inclination of the dipole to the observer, but for 
i = 45° the average field strength on the star is <~ 1.4 times 
the values given in the Table. Because of differences in 
underlying assumptions, these predictions are not identical, 
but they do have the same general dependence on system 
characteristics. Consequently, field strengths predicted by 
the 3 theories, while different in scale, nonetheless have the 
same pattern from star to star. Relatively weak fields are 
predicted for some stars (DN Tau, IP Tau), but detectably 
strong fields are expected on stars such as BP Tau. 



2.2. Measurement Techniques 

Virtually all measurements of stellar magnetic fields 
make use of the Zeeman effect. Typically, one of two gen- 
eral aspects of the Zeeman effect is utilized: (1) Zeeman 
broadening of magnetically sensitive lines observed in in- 
tensity spectra, or (2) circular polarization of magnetically 
sensitive lines. Due to the nature of the Zeeman effect, the 
splitting due to a magnetic field is proportional to A 2 of the 
transition. Compared with the A 1 dependence of Doppler 
line broadening mechanisms, this means that observations 
in the infrared (IR) are generally more sensitive to the pres- 
ence of magnetic fields than optical observations. 

The simplest model of the spectrum from a magnetic star 
assumes that the observed line profile can be expressed as 
F(X) = F B (X) * f + Fq(X) * (1 - /); where F B is the 
spectrum formed in magnetic regions, Fq is the spectrum 
formed in non-magnetic (quiet) regions, and / is the flux 
weighted surface filling factor of magnetic regions. The 
magnetic spectrum, Fb, differs from the spectrum in the 
quiet region not only due to Zeeman broadening of the line, 
but also because magnetic fields affect atmospheric struc- 
ture, causing changes in both line strength and continuum 
intensity at the surface. Most studies assume that the mag- 
netic atmosphere is in fact the same as the quiet atmosphere 
because there is no theory to predict the structure of the 
magnetic atmosphere. If the stellar magnetic field is very 
strong, the splitting of the a components is a substantial 
fraction of the line width, and it is easy to see the a compo- 
nents sticking out on either side of a magnetically sensitive 
line. In this case, it is relatively straightforward to measure 
the magnetic field strength, B. Differences in the atmo- 
spheres of the magnetic and quiet regions primarily affect 
the value of /. If the splitting is a small fraction of the in- 
trinsic line width, then the resulting observed profile is only 
subtly different from the profile produced by a star with no 
magnetic field and more complicated modelling is required 
to be sure all possible non-magnetic sources (e.g., rotation 
and pressure broadening) have been properly constrained. 

In cases where the Zeeman broadening is too subtle to 
detect directly, it is still possible to diagnose the presence of 
magnetic fields through their effect on the equivalent width 
of magnetically sensitive lines. For strong lines, the Zee- 
man effect moves the a components out of the partially sat- 
urated core into the line wings where they can effectively 
add opacity to the line and increase the equivalent width. 
The exact amount of equivalent width increase is a compli- 
cated function of the line strength and Zeeman splitting pat- 
tern (Basri et al, 1992). This method is primarily sensitive 
to the product of B multiplied by the filling factor / (Basri 
et al, 1992, Guenther et al, 1999). Since this method re- 
lies on relatively small changes in the line equivalent width, 
it is very important to be sure other atmospheric parameters 
which affect equivalent width (particularly temperature) are 
accurately measured. 

Measuring circular polarization in magnetically sensitive 
lines is perhaps the most direct means of detecting magnetic 
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Table 1 

Predicted Magnetic Field Strengths 
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AA Tau 


0.53 


1.74 


0.33 


8.20 


810 


240 


960 


8.0 


2.57 


BP Tau 


0.49 


1.99 


2.88 


7.60 


1370 


490 


1620 


6.4 


2.17 


CY Tau 


0.42 


1.63 


0.75 


7.90 


1170 


390 


1380 


7.7 




DE Tau 


0.26 


2.45 


2.64 


7.60 


420 


164 


490 


4.2 


1.35 


DF Tau 


0.27 


3.37 


17.7 


8.50 


490 


220 


570 


3.4 


2.98 


DK Tau 


0.43 


2.49 


3.79 


8.40 


810 


300 


950 


5.3 


2.58 


DN Tau 


0.38 


2.09 


0.35 


6.00 


250 


80 


300 


4.8 


2.14 


GG Tau A 


0.44 


2.31 


1.75 


10.30 


890 


320 


1050 


6.6 


1.57 


GI Tau 


0.67 


1.74 


0.96 


7.20 


1450 


450 


1700 


7.9 


2.69 


GK Tau 


0.46 


2.15 


0.64 


4.65 


270 


90 


320 


4.2 


2.13 


GM Aur 


0.52 


1.78 


0.96 


12.00 


1990 


660 


2340 


10.0 




IP Tau 


0.52 


1.44 


0.08 


3.25 


240 


60 


280 


5.2 




TWHya 


0.70 


1.00 


0.20 


2.20 


900 


240 


1060 


6.3 


2.61 


T Tau 


2.11 


3.31 


4.40 


2.80 


390 


110 


460 


3.2 


2.39 



Note. — Magnetic field values come from applying the theory of (a) Konigl (1991), (6) Cameron and 
Campbell, or (c) Shu et al. (1994). These are the equatorial field strengths assuming a dipole magnetic 
field. 



fields on stellar surfaces, but is also subject to several limi- 
tations. When viewed along the axis of a magnetic field, the 
Zeeman a components are circularly polarized, but with op- 
posite helicity; and the 7r component is absent. The helicity 
of the a components reverses as the polarity of the field re- 
verses. Thus, on a star like the Sun that typically displays 
equal amounts of + and — polarity fields on its surface, the 
net polarization is very small. If one magnetic polarity does 
dominate the visible surface of the star, net circular polar- 
ization is present in Zeeman sensitive lines, resulting in a 
wavelength shift between the line observed through right- 
and left-circular polarizers. The magnitude of the shift rep- 
resents the surface averaged line of sight component of the 
magnetic field (which on the Sun is typically less than 4 G 
even though individual magnetic elements on the solar sur- 
face range from ~ 1.5 kG in plages to ~ 3.0 kG in spots). 
Several polarimetric studies of cool stars have generally 
failed to detect circular polarization, placing limits on the 
disk-averaged magnetic field strength present of 10 — 100 
G (e.g., Vogt, 1980; Brown and Landstreet, 1981; Borra et 
al, 1984). One notable exception is the detection of circu- 
lar polarization in segments of the line profile observed on 
rapidly rotating dwarfs and RS CVn stars where Doppler 
broadening of the line "resolves" several independent strips 
on the stellar surface (e.g., Donati et al., 1997; Petit et al, 
2004; Jardine et al, 2002). 



2.3. Mean Magnetic Field Strength 

TTS typically have wsini values of 10 km s _1 , which 
means that observations in the optical typically cannot de- 
tect the actual Zeeman broadening of magnetically sensitive 
lines because the rotational broadening is too strong. Never- 
theless, optical observations can be used with the equivalent 
width technique to detect stellar fields. Basri et al. (1992) 
were the first to detect a magnetic field on the surface of a 
TTS, inferring a value of Bf = 1.0 kG on the NTTS Tap 
35. For the NTTS Tap 10, Basri et al. (1992) find only an 
upper limit of Bf < 0.7 kG. Guenther et al. (1999) apply 
the same technique to spectra of 5 TTSs, claiming signif- 
icant field detections on two stars; however, these authors 
analyze their data using models off by several hundred K 
from the expected effective temperature of their target stars, 
always a concern when relying on equivalent widths. 

As we saw above, observations in the IR will help solve 
the difficulty in detecting direct Zeeman broadening. For 
this reason and given the temperature of most TTSs (K7 - 
M2), Zeeman broadening measurements for these stars are 
best done using several Ti I lines found in the K band. Ro- 
bust Zeeman broadening measurements require Zeeman in- 
sensitive lines to constrain nonmagnetic broadening mech- 
anisms. Numerous CO lines at 2.31 ^m have negligible 
Lande-g factors, making them an ideal null reference. 

It has now been shown that the Zeeman insensitive CO 
lines are well fitted by models with the same level of rota- 
tional broadening as that determined from optical line pro- 
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files {Johns-Krull and Volenti, 2001; Johns-Krull et al, 
2004; Yang et al, 2005). In contrast, the 2.2 /im Ti I lines 
cannot be fitted by models without a magnetic field. In- 
stead, the observed spectrum is best fit by a model with 
a superposition of synthetic spectra representing different 
regions on the star with different magnetic field strengths. 
Typically, the field strengths in these regions are assumed to 
have values of 0, 2, 4, and 6 kG and only the filling factor of 
each region is solved for. The resulting magnetic field dis- 
tribution is unique because the Zeeman splitting produced 
by a 2 kG field is comparable to the nonmagnetic width of 
the Ti I spectral lines. In other words, the Zeeman resolu- 
tion of the Ti I lines is about 2 kG (see Johns-Krull et al, 
1999b, 2004). 

The intensity-weighted mean magnetic field strength, B, 
over the entire surface of most TTSs analyzed to date is 
~ 2.5 kG, with field strengths reaching at least 4 kG and 
probably even 6 kG in some regions. Thus, magnetic fields 
on TTSs are stronger than on the Sun, even though the 
surface gravity on these stars is lower by a factor of ten. 
On the Sun and other main-sequence stars, magnetic field 
strength seems to be set by an equipartition of gas and mag- 
netic pressure. In contrast, the photospheres of TTSs are 
apparently dominated by magnetic pressure, rather than gas 
pressure (see also Johns-Krull et al, 2004). Strong mag- 
netic fields are ubiquitous on TTSs. By fitting IR spectra, 
magnetic field distributions for several TTSs have now been 
measured (Johns-Krull et al., 1999b, 2001, 2004; Yang et 
al., 2005). Many of these field strengths are reported in Ta- 
ble 1. 

2.4. Magnetic Field Topology 

Zeeman broadening measurements are sensitive to the 
distribution of magnetic field strengths, but they have lim- 
ited sensitivity to magnetic geometry. In contrast, circular 
polarization measurements for individual spectral lines are 
sensitive to magnetic geometry, but they provide limited in- 
formation about field strength. The two techniques comple- 
ment each other well, as we demonstrate below. 

Most existing magnetospheric accretion models assume 
that intrinsic TTS magnetic fields are dipolar, but this would 
be unprecedented for cool stars. The higher order compo- 
nents of a realistic multi-polar field will fall off more rapidly 
with distance than the dipole component, so at the inner 
edge of the disk a few stellar radii from the surface, it is 
likely the dipolar component of the stellar field will dom- 
inate. However, at the stellar surface the magnetic field is 
likely to be more complicated. In support of this picture 
is the fact that spectropolarimetric observations do not de- 
tect polarization in photospheric absorption lines: Brown 
and Landstreet (1981) failed to detect polarization in T Tau 
and two FU Ori objects; Johnstone and Penston (1986) ob- 
served 3 CTTSs and reported a marginal field detection for 
RU Lup, but they were not able to confirm the signal in a 
subsequent observation, perhaps because of rotational mod- 
ulation (Johnst one and Penston, \9%1); Donati et al. (1997) 
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Fig. 2. — Variations in the circular polarization of the He I 
emission line as a function of rotation phase for 4 CTTSs. 
Polarization levels are translated into B z values in the line 
formation region. Vertical bars centered on each measure- 
ment (x) give the la uncertainty in the field measurement. 
Solid lines show predicted rotational modulation in B z for 
a single magnetic spot at latitudes ((f)) ranging from 0° to 
90° in 15° increments. The best fit latitude is shown in the 
thick solid curve. 

find no evidence for a strong dipolar field component in the 
3 TTSs they observed; Johns-Krull et al. (1999a) failed 
to detect polarization in the photospheric lines of BP Tau; 
and Valenti and Johns-Krull (2004) do not detect signifi- 
cant polarization in the photospheric lines of 4 CTTSs each 
observed over a rotation period. Smirnov et al. (2003) re- 
port a marginal detection of circular polarization in the lines 
of T Tau corresponding to a field of ~ 150 ± 50 G; how- 
ever, Smirnov et al. (2004) and Daou et al. (2005) failed to 
confirm this detection, placing an upper limit on the field of 

< 120 G for T Tau. 

However, Johns-Krull et al. (1999a) did discover cir- 
cular polarization in CTTS emission line diagnostics that 
form predominantly in the accretion shock at the surface of 
the star. This circular polarization signal is strongest in the 
narrow component of the He 1 5876 A emission line, but it is 
also present in the Ca I infrared triplet lines. The peak value 
of B z is 2.5 kG, which is comparable to our measured val- 
ues of B. Circular polarization in the He I 5876 A emission 
line has now been observed in a number of CTTSs (Valenti 
et al., 2004; Symington et al, 2005b). Note, since this po- 
larization is detected in a line associated with the accretion 
shock on CTTSs, it forms over an area covering typically 

< 5% of the stellar surface (Valenti et al., 1993; Calvet and 
Gullbring, 1998). While the field in this 5% of the star ap- 
pears to be highly organized (and as discussed below may 
trace the dipole component of the field at the surface), the 
lack of polarization detected in photospheric lines forming 
over the entire surface of the star strongly rule out a global 
dipole geometry for the entire field. 

Figure 2 shows measurements of B z on 6 consecutive 
nights. These measurements were obtained at McDon- 
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aid Observatory, using the Zeeman analyzer described by 
Johns-Krull et al. (1999a). The measured values of B z 
vary smoothly on rotational timescales, suggesting that uni- 
formly oriented magnetic field lines in accretion regions 
sweep out a cone in the sky, as the star rotates. Rotational 
modulation implies a lack of symmetry about the rotation 
axis in the accretion or the magnetic field or both. For ex- 
ample, the inner edge of the disk could have a concentra- 
tion of gas that corotates with the star, preferentially illu- 
minating one sector of a symmetric magnetosphere. Al- 
ternatively, a single large scale magnetic loop could draw 
material from just one sector of a symmetric disk. 

Figure 2 shows one interpretation of the He I polariza- 
tion data. Predicted values of B z are shown for a simple 
model consisting of a single magnetic spot at latitude <fi that 
rotates with the star. The magnetic field is assumed to be 
radial with a strength equal to our measured values of B. 
Inclination of the rotation axis is constrained by measured 
v sin i and rotation period, except that inclination (i) is al- 
lowed to float when it exceeds 60° because v sin i measure- 
ments cannot distinguish between these possibilities. Pre- 
dicted variations in B z are plotted for spot latitudes ranging 
from 0° to 90° in 15° increments. The best fitting model is 
shown by the thick curve. The corresponding spot latitude 
and reduced x 2 are given on the right side of each panel. 
The null hypothesis (that no polarization signal is present) 
produces very large values of \ 2 which are given on the left 
side of each panel. In all four cases, this simple magnetic 
spot model reproduces the observed B z time series. The 
He I rotationally modulated polarization combined with the 
lack of detectable polarization in photospheric absorption 
lines as described above paints a picture in which the mag- 
netic field on TTSs displays a complicated geometry at the 
surface which gives way to a more ordered, dipole-like ge- 
ometry a few stellar radii from the surface where the field 
intersects the disk. The complicated surface topology re- 
sults in no net polarization in photospheric absorption lines, 
but the dipole-like geometry of the field at the inner disk 
edge means that accreting material follows these field lines 
down to the surface so that emission lines formed in the 
accretion shock preferentially illuminate the dipole compo- 
nent of the field, producing substantial circular polarization 
in these emission lines. 

2.5. Confronting Theory with Observations 

At first glance, it might appear that magnetic field mea- 
surements on TTS are generally in good agreement with 
theoretical expectations. Indeed, the IR Zeeman broaden- 
ing measurements indicate mean fields on several TTSs of 
~ 2 kG, similar in value to those predicted in Table 1 (recall 
the field values in the Table are the equatorial values for a 
dipolar field, and that the mean field is about 1 .5 times these 
equatorial values). However, in detail the field observations 
do not agree with the theory. This can be seen in Figure 3, 
where we plot the measured magnetic field strengths ver- 
sus the predicted field strengths from Shu et al. (1994, see 
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Fig. 3. — Observed mean magnetic field strength deter- 
mined from IR Zeeman broadening measurements as a 
function of the predicted field strength from Table 1 for the 
theory of Shu et al. (1994). No statistically significant cor- 
relation is found between the observed and predicted field 
strengths. 

Table 1). Clearly, the measured field strengths show no cor- 
relation with the predicted field strengths. The field topol- 
ogy measurements give some indication to why there may 
be a lack of correlation: the magnetic field on TTSs are not 
dipolar, and the dipole component to the field is likely to 
be a factor of <~ 10 or more lower than the values predicted 
in Table 1. As discussed in Johns-Krull et al. (1999b), 
the 3 studies which produce the field predictions in Table 1 
involve uncertain constants which describe the efficiency 
with which stellar field lines couple to the accretion disk. 
If these factors are much different than estimated, it may be 
that the required dipole components to the field are substan- 
tially less than the values given in the Table. On the other 
hand, equation 2.2 was derived assuming perfect coupling 
of the field and the matter, so it serves as a firm upper limit 
to Rt as discussed in §2.1. Spectropolarimetry of TTSs in- 
dicates that the dipole component of the magnetic field is 
< 0.1 kG (Volenti and Johns-Krull, 2004; Smirnov et al., 
2004; Daou et al, 2005). Putting this value into equation 
2.2, we find Rt < 1.9 R* for typical CTTS parameters. 
Such a low value for the truncation radius is incompatible 
with rotation periods of 7- 1 days as found for many CTTSs 
(Table 1 and, e.g., Herbst et al, 2002). 

Does this then mean that magnetospheric accretion does 
not work? Independent of the coupling efficiency between 
the stellar field and the disk, magnetospheric accretion mod- 
els predict correlations between stellar and accretion param- 
eters. As shown in §2.3, the fields on TTS are found to all 
be rather uniform in strength. Eliminating the stellar field 
then, Johns-Krull and Gafford (2002) looked for correla- 
tion among the stellar and accretion parameters, finding lit- 
tle evidence for the predicted correlations. This absence of 
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the expected correlations had been noted earlier by Muze- 
rolle et al. (2001). On the other hand, Johns-Krull and 
Gafford (2002) showed how the models of Ostriker and Shu 
(1995) could be extended to take into account non-dipole 
field geometries. Once this is done, the current data do re- 
veal the predicted correlations, suggesting magnetospheric 
accretion theory is basically correct as currently formulated. 
So then, how do we reconcile the current field measure- 
ments with this picture? While the dipole component of 
the field is small on TTSs, it is clear the stars posses strong 
fields over most, if not all, of their surface but with a com- 
plicated surface topology. Perhaps this can lead to a strong 
enough field so that Rt ~ 6 i?* as generally suggested by 
observations of CTTS phenomena. More complicated nu- 
merical modelling of the interaction of a complex geometry 
field with an accretion disk will be required to see if this is 
feasible. 



3. SPECTRAL DIAGNOSTICS OF MAGNETO- 
SPHERIC ACCRETION 

Permitted emission line profiles from CTTSs, in partic- 
ular the Balmer series, show a wide variety of morpholo- 
gies including symmetric, double-peaked, P Cygni, and in- 
verse P Cygni (IPC) type {Edwards et al, 1994): com- 
mon to all shapes is a characteristic line width indicative of 
bulk motion within the circumstellar material of hundreds 
of kms -1 . The lines themselves encode both geometrical 
and physical information on the accretion process and its 
rate, and the challenge is to use the profiles to test and re- 
fine the magnetospheric accretion model. 

Interpretation of the profiles requires a translational step 
between the physical model and the observable spectra; this 
is the process of radiative-transfer (RT) modelling. The 
magnetospheric accretion paradigm presents a formidable 
problem in RT, since the geometry is two or three dimen- 
sional, the material is moving, and the radiation-field and 
the accreting gas are decoupled (i.e. the problem is non- 
LTE). However, the past decade has seen the development 
of increasingly sophisticated RT models that have been used 
to model line profiles (both equivalent width and shape) in 
order to determine accretion rates. In this section we de- 
scribe the development of these models, and characterize 
their successes and failures. 

Current models are based on idealized axisymmetric 
geometry, in which the circumstellar density structure is 
calculated assuming free-fall along dipolar field lines that 
emerge from a geometrically thin disc at a range of radii 
encompassing the corotation radius. It is assumed that the 
kinetic energy of the accreting material is completely ther- 
malized, and that the accretion luminosity, combined with 
the area of the accretion footprints (rings) on the stellar sur- 
face, provide the temperature of the hot spots. The circum- 
stellar density and velocity structure is then fully described 
by the mass accretion rate, and the outer and inner radii of 



the magnetosphere in terms of the photospheric radius of 
the star (Hartmann et al, 1994). 

A significant, but poorly constrained, input parameter for 
the models is the temperature structure of the accretion flow. 
This is a potential pitfall, as the form of the temperature 
structure may have a significant impact on the line source 
functions, and therefore the line profiles themselves. Self- 
consistent radiative equilibrium models {Martin, 1996) in- 
dicate that adiabatic heating and cooling via bremsstrahlung 
dominate the thermal budget, whereas Hartmann and co- 
workers adopt a simple volumetric heating rate combined 
with a schematic radiative cooling rate which leads to a 
temperature structure that goes as the reciprocal of the den- 
sity. Thus the temperature is low near the disc, and passes 
through a maximum (as the velocity increases and density 
decreases) before the stream cools again as it approaches 
the stellar surface (and the density increases once more). 

With the density, temperature and velocity structure of 
the accreting material in place, the level populations of the 
particular atom under consideration must be calculated un- 
der the constraint of statistical equilibrium. This calcula- 
tion is usually performed using the Sobolev approximation, 
in which it is assumed that the conditions in the gas do not 
vary significantly over a length scale given by 

IS = Vtherm/(dv/dr) (1) 

where v t herm is the thermal velocity of the gas and dv/dr is 
the velocity gradient. Such an approximation is only strictly 
valid in the fastest parts of the accretion flow. Once the level 
populations have converged, the line opacities and emissiv- 
ities are then computed, allowing the line profile of any par- 
ticular transition to be calculated. 

The first models computed using the method outlined 
above were presented by Hartmann et al (1994), who 
adopted a two-level atom approximation. It was demon- 
strated that the magnetospheric accretion model could re- 
produce the main characteristics of the profiles, including 
IPC profiles and blue-shifted central emission peaks. The 
original Hartmann et al. model was further improved by 
Muzerolle et al. (2001). Instead of using a two-level ap- 
proximation, they solved statistical equilibrium (still un- 
der Sobolev) for a 20-level hydrogen atom. Their line 
profiles were computed using a direct integration method, 
which, unlike the Sobolev approach, allows the inclusion 
of Stark broadening effects. It was found that the broad- 
ening was most significant for Ha, with the line reaching 
widths of <~ 500 kms -1 , a width that significantly exceeds 
the doppler broadening due to infall alone, and in much 
better agreement with observation. The H/3 model profiles 
were found to be in broad agreement with the observations, 
in terms of the velocity of the emission peak (Alencar and 
Basri, 2000), and in the asymmetry of the profiles (Edwards 
et al, 1994). Figure 4 shows model Ha profiles as a func- 
tion of mass accretion rate and accretion flow temperature; 
one can see that for typical CTTS accretion rates the line 
profiles are broadly symmetric although slightly blueshifted 
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Fig. 4. — Ha model profiles for a wide range of mass ac- 
cretion rate and accretion flow maximum temperature (from 
Kurosawa et ai, 2006). The profiles are based on canonical 
CTTS parameters (R = 2 R Q , M = 0.5 M , T = 4000 K) 
viewed at an inclination of 55°. The maximum temperature 
of the accretion flow is indicated along the left of the figure, 
while the accretion rate (in M Q yr _1 ) is shown along the 
top. 

- the reduced optical depth for the lower accretion rate mod- 
els yields the IPC morphology. 

Axisymmetric models are obviously incapable of repro- 
ducing the wide range of variability that is observed in the 
emission lines of CTTSs (Sect. 4). Although the addition 
of further free parameters to models naturally renders them 
more arbitrary, the observational evidence for introducing 
such parameters is compelling. Perhaps the simplest ex- 
tension is to break the axisymmetry of the dipole, leaving 
curtains of accretion in azimuth - models such as these have 
been proposed by a number of observers attempting to ex- 
plain variability in CTTSs and are observed in MHD simu- 
lations (Romanova et al, 2003). Synthetic time-series for a 
CTTS magnetosphere structured along these lines were pre- 
sented by Symington et al. (2005a). It was found that some 
gross characteristics of the observed line profiles were pro- 
duced using a 'curtains' model, although the general level 
of variability predicted is larger than that observed, suggest- 
ing that the magnetosphere may be characterized by a high 
degree of axisymmetry, broken by higher-density streams 
that produce the variability. 

The emission line profiles of CTTSs often display the 
signatures of outflow as well as infall, and recent attempts 
have been made to account for this in RT modelling. Alen- 
car et al. (2005) investigated a dipolar accretion geome- 
try combined with a disk wind in order to model the line 
profile variability of RW Aur. They discovered that magne- 
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Fig. 5. — Sample Ha model profiles (Kurosawa et al., 
2006) which characterize the morphological classification 
(Types I-IV B/R) by Reipurth et al. (1996). The combina- 
tion of magnetospheric accretion, the accretion disc, and the 
collimated disk wind can reproduce the wide range of Ha 
profiles seen in observations. The horizontal axes are ve- 
locities in kms -1 and the vertical axes are continuum nor- 
malized intensities. 

tospheric accretion alone could not simultaneously model 
Ha, H/3 and NaD profiles, and found that the wind contri- 
bution to the lines profiles is quite important in that case. 

Hybrid models (Kurosawa et al, 2006) combining a 
standard dipolar accretion flow with an outflow (e.g., Fig. 6) 
are capable of reproducing the broad range of observed Ha 
profiles (Fig. 5). Obviously spectroscopy alone is insuffi- 
cient to uniquely identify a set of model parameters for an 
individual object, although by combining spectroscopy with 
other probes of the circumstellar material, one should be 
able to reduce the allowable parameter space considerably. 
For example linear spectropolarimetry provides a unique in- 
sight into the accretion process; scattering of the line emis- 
sion by circumstellar dust imprints a polarization signature 
on the line which is geometry dependent. An Ha spectropo- 
larimetric survey by Vink et al. (2005a) revealed that 9 out 
of 1 1 CTTSs showed a measurable change in polarization 
through the line, while simple numerical models by Vink et 
al. (2005b) demonstrate that this polarization may be used 
to gauge the size of the disk inner hole. 

The radiative-transfer models described above are now 
routinely used to determine mass accretion rates across the 
mass spectrum from Herbig AeBe (Muzerolle et al, 2004) 
stars to brown dwarfs (Lawson et al, 2004; Muzerolle et 
al, 2005), and in the CTTS mass regime at least the ac- 
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Fig. 6. — A simulated Ha image of an accreting CTTS with 
an outflow (logM acc = —8, logM W j n d = —9) viewed 
at an inclination of 80°. The wind emission is negligible 
compared to the emission from the magnetosphere, and the 
lower half of the wind is obscured by the circumstellar disk 
(Kurosawa et at, 2006). 



cretion rates derived from RT modelling have been roughly 
calibrated against other accretion-rate measures, such as the 
UV continuum (e.g., Muzerolle et al, 2001). However, one 
must be aware of the simplifying assumptions which un- 
derlie the models and that must necessarily impact on the 
validity of any quantity derived from them, particularly the 
mass accretion rate. Magnetic field measurements (Sect. 2) 
and time-series spectroscopy (Sect. 4) clearly show us that 
the geometry of the magnetosphere is far from a pristine 
axisymmetric dipole, but instead probably consists of many 
azimuthally distributed funnels of accretion, curved by ro- 
tation and varying in position relative to the stellar surface 
on the timescale of a few stellar rotation periods. Further- 
more, the temperature of the magnetosphere and the mass 
accretion rate are degenerate quantities in the models, with 
a higher temperature magnetosphere producing more line 
flux for the same accretion rate. This means that brown 
dwarf models require a much higher accretion stream tem- 
perature than those of CTTSs in order to produce the ob- 
served line flux, and although the temperature is grossly 
constrained by the line broadening (which may preclude 
lower temperature streams) the thermal structure of the ac- 
cretion streams is still a problem. Despite these uncertain- 
ties, and in defense of the BD models, it should be noted 
that the low accretion rates derived are consistent with both 
the lack of optical veiling (Muzerolle et al, 2003a) and the 
strength of the Ca II A8662 line (Mohanty et al, 2005). 

Current models do not match the line core particu- 
larly well, which is often attributed to a break down of 
the Sobolev approximation; co-moving frame calculations 
(which are many orders of magnitude more expensive com- 
putationally) may be required. An additional problem with 



current RT modelling is the reliance on fitting a single pro- 
file - current studies have almost always been limited to 
Ha - one that rarely shows an IPC profile (Edwards et al, 
1994; Reipurth et al, 1996), is vulnerable to contamination 
by outflows (e.g., Alencar et al, 2005) and may be signif- 
icantly spatially extended (Takami et al, 2003). Even in 
modelling a single line, it is fair to say that the state-of-the- 
art is some way short of line profile fitting; the best fits re- 
ported in the literature may match the observation in terms 
of peak intensity, equivalent width, or in the line wings, 
but are rarely convincing reproductions of the observations 
in detail. Only by simultaneously fitting several lines may 
one have confidence in the models, particularly if those 
lines share a common upper/lower level (Ha and Pa/3 for 
example). Although such observations are in the literature 
(e.g., Edwards et al, 1994; Folha and Emerson, 2001) their 
usefulness is marginalized by the likely presence of signif- 
icant variability between the epochs of the observations at 
the different wavelengths: simultaneous observations of a 
wide range of spectral diagnostics are required. Despite 
the caveats described above, line profile modelling remains 
a useful (and in the BD case the only) route to the mass 
accretion rate, and there is real hope that the current fac- 
tor ~ 5 uncertainties in mass accretion rates derived from 
RT modelling may be significantly reduced in the future. 



4. OBSERVATIONAL EVIDENCE FOR MAGNE- 
TOSPHERIC ACCRETION 

Observations seem to globally support the magneto- 
spheric accretion concept in CTTSs, which includes the 
presence of strong stellar magnetic fields, the existence of 
an inner magnetospheric cavity of a few stellar radii, mag- 
netic accretion columns filled with free falling plasma, and 
accretion shocks at the surface of the stars. While this sec- 
tion summarizes the observational signatures of magneto- 
spheric accretion in T Tauri stars, there is some evidence 
that the general picture applies to a much wider range of 
mass, from young brown dwarfs (Muzerolle et al, 2005; 
Mohanty et al, 2005) to Herbig Ae-Be stars (Muzerolle et 
al, 2004; Calvet et al, 2004; Sorelli et al, 1996). 

In recent years, the rapidly growing number of detec- 
tions of strong stellar magnetic fields at the surface of young 
stars seem to put the magnetospheric accretion scenario on a 
robust ground (see Sect. 2). As expected from the models, 
given the typical mass accretion rates (10~ 9 to 10~ 7 M Q 
yr _1 , Gullbring et al, 1998) and magnetic field strengths 
(2 to 3 kG, Volenti and Johns-Krull, 2004) obtained from 
the observations, circumstellar disk inner holes of about 3- 
9 R* are required to explain the observed line widths of 
the CO fundamental emission, that likely come from gas in 
Keplerian rotation in the circumstellar disk of CTTSs (Na- 
jita et al, 2003). There has also been evidence for accre- 
tion columns through the common occurrence of inverse P 
Cygni profiles with redshifted absorptions reaching several 
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hundred km s 1 , which indicates that gas is accreted onto 
the star from a distance of a few stellar radii (Edwards et 
al, 1994). 

Accretion shocks are inferred from the rotational mod- 
ulation of light curves by bright surface spots (Bouvier et 
al, 1995) and modelling of the light curves suggests hot 
spots covering about one percent of the stellar surface. The 
theoretical prediction of accretion shocks and its associated 
hot excess emission are also supported by accretion shock 
models that successfully reproduce the observed spectral 
energy distributions of optical and UV excesses (Calvet and 
Gullbring, 1998; Ardila and Basri, 2000; Gullbring et al., 
2000). In these models, the spectral energy distribution of 
the excess emission is explained as a combination of opti- 
cally thick emission from the heated photosphere below the 
shock and optically thin emission from the preshock and 
postshock regions. Gullbring et al. (2000) also showed 
that the high mass accretion rate CTTSs have accretion 
columns with similar values of energy flux as the moder- 
ate to low mass accretion rate CTTSs, but their accretion 
columns cover a larger fraction of the stellar surface (filling 
factors ranging from less than 1 % for low accretors to more 
than 10% for the high one). A similar trend was observed 
by Ardila and Basri (2000) who found, from the study of 
the variability of IUE spectra of BP Tau, that the higher the 
mass accretion rate, the bigger the hot spot size. 

Statistical correlations between line fluxes and mass ac- 
cretion rates predicted by magnetospheric accretion mod- 
els have also been reported for emission lines in a broad 
spectral range, from the UV to the near-IR {Johns-Krull et 
al., 2000; Beristain et al., 2001; Alencar and Basri, 2000; 
Muzerolle et al, 2001; Folha and Emerson, 2001). How- 
ever, in recent years, a number of observational results in- 
dicate that the idealized steady-state axisymmetric dipolar 
magnetospheric accretion models cannot account for many 
observed characteristics of CTTSs. 

Recent studies showed that accreting systems present 
strikingly large veiling variability in the near-IR (Eiroa et 
al, 2002; Barsony et al, 2005), pointing to observational 
evidence for time variable accretion in the inner disk. More- 
over, the near-IR veiling measured in CTTSs is often larger 
than predicted by standard disk models (Folha and Emer- 
son, 1999; Johns-Krull and Volenti, 2001). This suggests 
that the inner disk structure is significantly modified by its 
interaction with an inclined stellar magnetosphere and thus 
departs from a flat disk geometry. Alternatively, a "puffed" 
inner disk rim could result from the irradiation of the in- 
ner disk by the central star and accretion shock (Natta et al. 
2001; Muzerolle et al, 2003b). In mildly accreting T Tauri 
stars, the dust sublimation radius computed from irradiation 
models is predicted to lie close to the corotation radius (3-9 
R* ~ 0.03-0.08 AU) though direct interferometric measure- 
ments tend to indicate larger values (0.08-0.2 AU, Akeson 
etal, 2005). 

Observational evidence for an inner disk warp has been 
reported by Bouvier et al. (1999, 2003) for AA Tau, as 
expected from the interaction between the disk and an in- 
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Fig. 7. — The rotational modulation of the Ha line profile 
of the CTTS AA Tau (8. 2d period). Line profiles are or- 
dered by increasing rotational phase (top panel number) at 
different Julian dates (bottom panel number). Note the de- 
velopment of a high velocity redshifted absorption compo- 
nent in the profile from phase 0.39 to 0.52, when the funnel 
flow is seen against the hot accretion shock (from Bouvier 
et al, in prep.). 

dined stellar magnetosphere (see Sect. 5). Inclined mag- 
netospheres are also necessary to explain the observed pe- 
riodic variations over a rotational timescale in the emis- 
sion line and veiling fluxes of a few CTTSs (Johns and 
Basri, 1995; Petrov et al, 1996; 2001; Bouvier et al, 1999; 
Batalha et al, 2002). These are expected to arise from the 
variations of the projected funnel and shock geometry as the 
star rotates. An example can be seen in Fig. 7 that shows the 
periodic modulation of the Ha line profile of the CTTS AA 
Tau as the system rotates, with the development of a high 
velocity redshifted absorption component when the funnel 
flow is seen against the hot accretion shock. Sometimes, 
however, multiple periods are observed in the line flux vari- 
ability and their relationship to stellar rotation is not al- 
ways clear (e.g., Alencar and Batalha, 2002; Oliveira et al, 
2000). The expected correlation between the line flux from 
the accretion columns, and the continuum excess flux from 
the accretion shock is not always present either (Ardila and 
Basri, 2000; Batalha et al, 2002), and the correlations pre- 
dicted by static dipolar magnetospheric accretion models 
are generally not seen (Johns-Krull and Gafford, 2002). 

Winds are generally expected to be seen as forbidden 
emission lines or the blueshifted absorption components of 
permitted emission lines. Some permitted emission line 
profiles of high-mass accretion rate CTTSs, however, do not 
always look like the ones calculated with magnetospheric 
accretion models and this could be in part due to a strong 
wind contribution to the emission profiles, given the high 
optical depth of the wind in these cases (Muzerolle et al, 
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2001; Alencar et al, 2005). Accretion powered hot winds 
originating at or close to the stellar surface have recently 
been proposed to exist in CTTSs with high mass accre- 
tion rates {Edwards et al, 2003). These winds are inferred 
from the observations of P Cygni profiles of the He I line 
(10780 A) that present blueshifted absorptions which ex- 
tend up to -400 km/s. Matt and Pudritz (2005) have argued 
that such stellar winds can extract a significant amount of 
the star's angular momentum, thus helping regulate the spin 
of CTTSs. Turbulence could also be important and help ex- 
plain the very wide (± 500 km s _1 ) emission line profiles 
commonly observed in Balmer and Mgll UV lines (Ardila 
etal, 2002). 

Synoptic studies of different CTTSs highlighted the dy- 
namical aspect of the accretion/ejection processes, which 
only recently has begun to be studied theoretically by nu- 
merical simulations (see Sect. 5). The accretion process 
appears to be time dependent on several timescales, from 
hours for non-steady accretion (Gullbring et al., 1996; 
Alencar and Batalha, 2002; Stempels and Piskunov, 2002; 
Bouvier et al, 2003) to weeks for rotational modulation 
{Smith et al, 1999; Johns and Basri, 1995; Petrov et al, 
2001), and from months for global instabilities of the mag- 
netospheric structure {Bouvier et al, 2003) to years for 
EXor and FUor eruptions (e.g., Reipurth and Aspin, 2004; 
Herbig, 1989). 

One reason for such a variability could come from the 
interaction between the stellar magnetosphere and the inner 
accretion disk. In general, magnetospheric accretion mod- 
els assume that the circumstellar disk is truncated close to 
the corotation radius and that field lines threading the disk 
corotate with the star. However, many field lines should in- 
teract with the disk in regions where the star and the disk 
rotate differentially. Possible evidence has been reported 
for differential rotation between the star and the inner disk 
(Oliveira et al, 2000) through the presence of an observed 
time delay of a few hours between the appearance of high 
velocity redshifted absorption components in line profiles 
formed in different regions of the accretion columns. This 
was interpreted as resulting from the crossing of an az- 
imuthally twisted accretion column on the line of sight. An- 
other possible evidence for twisted magnetic field lines by 
differential rotation leading to reconnection events has been 
proposed by Montmerle et al. (2000) for the embedded 
protostellar source YLW 15, based on the observations of 
quasi-periodic X-ray flaring. A third possible evidence was 
reported by Bouvier et al. (2003) for the CTTS AA Tau. On 
timescales of the order of a month, they observed significant 
variations in the line and continuum excess flux, indicative 
of a smoothly varying mass accretion rate onto the star. At 
the same time, they found a tight correlation between the 
radial velocity of the blueshifted (outflow) and redshifted 
(inflow) absorption components in the Ha emission line 
profile. This correlation provides support for a physical 
connection between time dependent inflow and outflow in 
CTTSs. Bouvier et al. (2003) interpreted the flux and ra- 
dial velocity variations in the framework of magnetospheric 



inflation cycles due to differential rotation between the star 
and the inner disk, as observed in recent numerical simula- 
tions (see Sect. 5). The periodicity of such instabilities, as 
predicted by numerical models, is yet to be tested observa- 
tionally and will require monitoring campaigns of chosen 
CTTSs lasting for several months. 

5. NUMERICAL SIMULATIONS OF MAGNETO- 
SPHERIC ACCRETION 

Significant progress has been made in recent years in the 
numerical modeling of magnetospheric accretion onto a ro- 
tating star with a dipolar magnetic field. One of the main 
problems is to find adequate initial conditions which do not 
destroy the disk in first few rotations of the star and do not 
influence the simulations thereafter. In particular, one must 
deal with the initial discontinuity of the magnetic field be- 
tween the disk and the corona, which usually leads to signif- 
icant magnetic braking of the disk matter and artificially fast 
accretion onto the star on a dynamical time-scale. Specific 
quasi-equilibrium initial conditions were developed, which 
helped to overcome this difficulty (Romanova et al, 2002). 
In axisymmetric (2D) simulations, the matter of the disk 
accretes inward slowly, on a viscous time-scale as expected 
in actual stellar disks. The rate of accretion is regulated 
by a viscous torque incorporated into the numerical code 
through the a prescription, with typically a v = 0.01 — 0.03. 

Simulations have shown that the accretion disk is dis- 
rupted by the stellar magnetosphere at the magnetospheric 
or truncation radius Rt, where the gas pressure in the disk 
is comparable to the magnetic pressure, P ram = B 2 /8tt 
(see Sect. 2). In this region matter is lifted above the disk 
plane due to the pressure force and falls onto the stellar 
surface supersonically along the field lines, forming fun- 
nel flows (Romanova et al, 2002). The location of the in- 
ner disk radius oscillates as a result of accumulation and 
reconnection of the magnetic flux at this boundary, which 
blocks or "permits" accretion (see discussion of this issue 
below), thus leading to non-steady accretion through the 
funnel flows. Nevertheless, simulations have shown that 
the funnel flow is a quasi-stationary feature during at least 
50 — 80 rotation periods of the disk at the truncation ra- 
dius, P , and recent simulations with improved numerical 
schemes indicate that this structure survives for more than 
1,000 Pq (Long et al, 2005). Axisymmetric simulations 
thus confirmed the theoretical ideas regarding the structure 
of the accretion flow around magnetized CTTSs. As a next 
step, similar initial conditions were applied to full 3D simu- 
lations of disk accretion onto a star with an inclined dipole, 
a challenging problem which required the development of 
new numerical methods (e.g., the "inflated cube" grid, cf. 
Koldoba et al, 2002; Romanova et al, 2003, 2004a). Sim- 
ulations have shown that the disk is disrupted at the trunca- 
tion radius Rt, as in the axisymmetric case, but the mag- 
netospheric flow to the star is more complex. Matter flows 
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Fig. 8. — A slice of the funnel stream obtained in 3D simulations for an inclined dipole (9 = 15°). The contour lines show 
density levels, from the minimum (dark) to the maximum (light). The corona above the disk has a low-density but is not 
shown. The thick lines depict magnetic field lines (from Romanova et al., 2004a). 



around the magnetosphere and falls onto the stellar surface 
supersonically. The magnetospheric structure varies de- 
pending on the misalignment angle of the dipole, but settles 
into a quasi-stationary state after a few P , as demonstrated 
by recent simulations run up to 40 Po (Kulkarni and Ro- 
manova, 2005). In both, 2D and 3D simulations the fluxes 
of matter and angular momentum to or from the star vary 
in time, however they are smooth on average. This average 
value is determined by the properties of the accretion disk. 

Numerical simulation studies have shown that a star may 
either spin up, spin down or be in rotational equilibrium 
when the net torque on the star vanishes. Detailed inves- 
tigation of the rotational equilibrium state has shown that 
the rotation of the star is then locked at an angular veloc- 
ity Q, eq which is smaller by a factor of ~ 0.67 — 0.83 than 
the angular velocity at the truncation radius {Long et ah, 
2005). The corresponding "equilibrium" corotation radius 
Rco ~ (1-3 — 1.5) Rt is close to that predicted theoret- 
ically (e.g., Ghosh and Lamb, 1978,1979b; Konigl, 1991). 
Recently, the disk-locking paradigm was challenged by a 
number of authors (e.g., Agapitou and Papaloizou, 2000; 
Matt and Pudritz, 2004, 2005). The skepticism was based 
on the fact that the magnetic field lines connecting the star 
to the disk may inflate and open, (e.g., Aly and Kuijpers, 
1990; Lovelace et al, 1995; Bardou, 1999; Uzdensky et al, 
2002), resulting in a significant decrease of angular momen- 
tum transport between the star and the disk. Such an open- 
ing of field lines was observed in a number of simulations 
(e.g., Miller and Stone, 1997; Romanova et al., 1998; Fendt 
and Elstner, 2000). Several factors, however, tend to restore 
an efficient disk-star connection. One of them is that the 
inflated field lines have a tendency to reconnect and close 
again (Uzdensky et al, 2002). Furthermore, there is always 



a region of closed field lines connecting the inner regions 
of the disk with the magnetosphere, which provides angu- 
lar momentum transport between the disk and the star (e.g., 
Pringle and Rees, 1972; Ghosh and Lamb, 1979b). This is 
the region where matter accretes through funnel flows and 
efficiently transports angular momentum to or from the star. 
This torque tends to bring a star in co-rotation with the in- 
ner regions of the disk. There is always, however, a smaller 
but noticeable negative torque either connected with the re- 
gion r > R co (Ghosh and Lamb, 1978,1979b), if the field 
lines are closed in this region, or associated to a wind which 
carries angular momentum out along the open field lines 
connecting the star to a low-density corona. Simulations 
have shown that the wind is magnetically-dominated (Long 
et al., 2005; Romanova et al., 2005), though the possibility 
of an accretion-driven stellar wind has also been discussed 
(Matt and Pudritz, 2005). The spin-down through magnetic 
winds was proposed earlier by Tout and Pringle (1992). 
Both torques are negative so that in rotational equilibrium 
a star rotates slower than the inner disk. Thus, the result 
is similar to the one predicted earlier theoretically, though 
the physics of the spin-down contribution may be different. 
Axisymmetric simulations of the fast rotating CTTSs have 
shown that they efficiently spin-down through both disk- 
magnetosphere interaction and magnetic winds (Romanova 
et al., 2005; Ustyugova et al, 2006). For instance, it was 
shown that a CTTS with an initial period P = 1 d spins 
down to the typically observed periods of about a week in 
less that 10 6 yr. 

Three-dimensional (3D) simulations of disk accretion 
onto a star with a misaligned dipolar magnetic field have 
shown that at the non-zero misalignment angle 9, where 9 
is an angle between the magnetic moment /i* and the rota- 
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= 30' 




Fig. 9. — 3D simulations show that matter accretes onto the 
star through narrow, high density streams (right panel) sur- 
rounded by lower density funnel flows that blanket nearly 
the whole magnetosphere (left panel). 




Fig. 10. — Top panels: matter flow close to the star at dif- 
ferent misalignment angles O. Bottom panels: the shape of 
the corresponding hot spots. Darker regions correspond to 
larger density. (From Romanova et al, 2004a). 

tional axis f2» of the star (with the disk axis aligned with 
n*), matter typically accretes in two and, in some cases, 
in several streams (Koldoba et al, 2002; Romanova et ai, 
2003, 2004a). Figure 8 shows a slice of the magnetospheric 
stream at = 15°. The density and pressure of the flow 
increase towards the star as a result of the convergence of 
the flow. They are also larger in the central regions of the 
funnel streams. Thus, the structure of the magnetospheric 
flow depends on the density. The high density part is chan- 
neled in narrow funnel streams, while the low density part is 
wider, with accreting matter blanketting the magnetosphere 
nearly completely {Romanova et ai, 2003, see Figure 9). 
The spectral lines which form in the funnel streams are red- 
shifted or blueshifted depending on the angle and view- 
ing angle i and their strength is modulated by the rotation 
of the star. 

Matter in the funnel flows falls onto the star's surface and 
forms hot spots. The shape of the spots and the distribution 
of different parameters (density, velocity, pressure) in the 
spots reflect those in the cross-section of the funnel streams 
(Romanova et ah, 2004a). Figure 10 shows an example of 
magnetospheric flows and hot spots at different 0. At rela- 
tively small angles, < 30°, the spots have the shape of a 
bow, while at very large angles, > 60°, they have a shape 



of a bar crossing the surface of the star near the magnetic 
pole. The density, velocity and pressure are the largest in 
the central regions of the spots and decrease outward (see 
Figure 10). The temperature also increases towards the cen- 
ter of the spots because the kinetic energy flux is the largest 
there. The rotation of the star with surface hot spots leads to 
variability with one or two peaks per period depending on 
and i. The two peaks are typical for larger and i. The po- 
sition of the funnel streams on the star is determined by both 
the angular velocity of the star and that of the inner radius of 
the disk. In the rotational equilibrium state, the funnel flows 
usually settle in a particular "favorite" position. However, if 
the accretion rate changes slightly, say, increases, then the 
truncation radius decreases accordingly and the angular ve- 
locity at the foot-point of the funnel stream on the disk is 
larger. As a result, the other end of the stream at the surface 
of the star changes its position by a small amount. Thus, 
the location of the spots "wobbles" around an equilibrium 
position depending on the accretion rate (Romanova et al., 
2004a). The variation of the accretion rate also changes the 
size and the brightness of the spots. 

The disk-magnetosphere interaction leads to the thicken- 
ing of the inner regions of the disk which eases the lifting of 
matter to the funnel flow. Matter typically accumulates near 
the closed magnetosphere forming a denser ring (Romanova 
et ai, 2002) which brakes into a spiral structure in case of 
misaligned dipole (Romanova et ai, 2003, 2004a). Typi- 
cally, two trailing spiral arms are obtained (see Figure 10). 
3D simulations have also shown that when accretion occurs 
onto a tilted dipole, the inner regions of the disk are slightly 
warped. This results from the tendency of disk material to 
flow along the magnetic equator of the misaligned dipole 
(Romanova et ai, 2003). Such a warping is observed for 
medium misalignment angles, 30° < < 60°. Disk warp- 
ing in the opposite direction (towards magnetic axis of the 
dipole) was predicted theoretically when the disk is strongly 
diamagnetic (Aly, 1980; Lipunov and Shakura, 1980; Lai, 
1999). The warping of the inner disk and the formation of a 
spiral structure in the accretion flow may possibly be at the 
origin of the observed variability of some CTTSs (Terquem 
and Papaloizou, 2000; Bouvier et ai, 2003). 

Progress has also been made in the modeling of outflows 
from the vicinity of the magnetized stars. Such outflows 
may occur from the disk-magnetosphere boundary (Shu et 
ai, 1994), from the disk (Blandford and Payne, 1982; Pu- 
dritz and Norman, 1986; Lovelace et ai, 1991; Lovelace 
et al., 1995; Casse and Ferreira, 2000; Pudritz et al., 
2006), or from the star (Matt and Pudritz, 2005). Magneto- 
centrifugally driven outflows were first investigated in pi- 
oneering short-term simulations by Hayashi et al. (1996) 
and Miller and Stone (1997) and later in longer-term sim- 
ulations with a fixed disk (Ouyed and Pudritz, 1997; Ro- 
manova et al., 1997; Ustyugova et al., 1999; Krasnopolsky 
et al., 1999; Fendt and Elstner, 2000). Simulations includ- 
ing feedback on the inner disk have shown that the process 
of the disk-magnetosphere interaction is non-stationary: the 
inner radius of the disk oscillates, and matter accretes to the 
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star and outflows quasi-periodically (Goodson et al, 1997, 
1999; Hirose et al, 1997; Matt et al, 2002; Kato et al, 
2004; Romanova et al, 2004b; Von Rekowski and Branden- 
burg, 2004; Romanova et al, 2005), as predicted by Aly 
and Kuijpers (1990). The characteristic timescale of vari- 
ability is determined by a number of factors, including the 
time-scale of diffusive penetration of the inner disk matter 
through the external regions of the magnetosphere {Good- 
son and Winglee, 1999). It was earlier suggested that re- 
connection of the magnetic flux at the disk-magnetosphere 
boundary may lead to X-ray flares in CTTSs (Hayashi et 
al, 1996; Feigelson and Montmerle, 1999) and evidence 
for very large flaring structures has been recently reported 
by Favata et al. (2005). 

So far simulations were done for a dipolar magnetic 
field. Observations suggest a non-dipolar magnetic field 
near the stellar surface (see Sect. 2, and also, e.g., Safier, 
1998; Kravtsova and Lamzin, 2003; Lamzin, 2003; Smirnov 
et al, 2005). If the dipole component dominates on the 
large scale, many properties of magnetospheric accretion 
will be similar to those described above, including the struc- 
ture of the funnel streams and their physical properties. 
However, the multipolar component will probably control 
the flow near the stellar surface, possibly affecting the shape 
and the number of hot spots. Simulations of accretion to a 
star with a multipolar magnetic field are more complicated, 
and should be done in the future. 

6. CONCLUSIONS 

Recent magnetic field measurements in T Tauri stars 
support the view that the accretion flow from the inner disk 
onto the star is magnetically controlled. While typical val- 
ues of 2.5 kG are obtained for photospheric fields, it also 
appears that the field topology is likely complex on the 
small-scales (R<R*), while on the larger scale (R>>R+) 
a globally more organized but weaker (<~0. 1 kG) magnetic 
component dominates. This structure is thought to inter- 
act with the inner disk to yield magnetically-channeled ac- 
cretion onto the star. Observational evidence for magneto- 
spheric accretion in classical T Tauri star is robust (inner 
disk truncation, hot spots, line profiles) and the rotational 
modulation of accretion/ejection diagnostics observed in 
some systems suggests that the stellar magnetosphere is 
moderately inclined relative to the star's rotational axis. Re- 
alistic 3D numerical models have capitalized on the obser- 
vational evidence to demonstrate that many properties of 
accreting T Tauri stars could be interpreted in the frame- 
work of magnetically-controlled accretion. One of the most 
conspicuous properties of young stars is their extreme vari- 
ability on timescales ranging from hours to months, which 
can sometimes be traced to instabilities or quasi-periodic 
phenomena associated to the magnetic star-disk interaction. 

Much work remains to be done, however, before reach- 
ing a complete understanding of this highly dynamical and 



time variable process. Numerical simulations still have to 
incorporate field geometries more complex than a tilted 
dipole, e.g., the superposition of a large-scale dipolar or 
quadrupolar field with multipolar fields at smaller scales. 
The modeling of emission line profiles now starts to com- 
bine radiative transfer computations in both accretion fun- 
nel flows and associated mass loss flows (disk winds, stellar 
winds), which indeed appears necessary to account for the 
large variety of line profiles exhibited by CTTSs. These 
models also have to address the strong line profile vari- 
ability which occurs on a timescale ranging from hours to 
weeks in accreting T Tauri stars. These foreseen devel- 
opments must be driven by intense monitoring of typical 
CTTSs on all timescales from hours to years, which com- 
bines photometry, spectroscopy and polarimetry in various 
wavelength domains. This will provide strong constraints 
on the origin of the variability of the various components of 
the star-disk interaction process (e.g., inner disk in the near- 
IR, funnel flows in emission lines, hot spots in the optical 
or UV, magnetic reconnections in X-rays, etc.). 

The implications of the dynamical nature of magneto- 
spheric accretion in CTTSs are plentiful and remain to be 
fully explored. They range from the evolution of stellar 
angular momentum during the pre-main sequence phase 
(e.g., Agapitou and Papaloizou, 2000), the origin of in- 
flow/ouflow short term variability (e.g., Woitas et ai, 2002; 
Lopez-Martin et al., 2003), the modeling of the near in- 
frared veiling of CTTSs and of its variations, both of which 
will be affected by a non planar and time variable inner disk 
structure (e.g., Carpenter et ai, 2001; Eiroa et al, 2002), 
and possibly the halting of planet migration close to the star 
(Linetal, 1996). 
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